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//for each pixel in the NxM image 
for(x=l: x<N; x+=l) 

{ 

for(y=l;y<M:y+=l) 
{ 

// Consider I(x,y), I(x-l,y), I(x,y-1),... neighborhood of (x,y) 

//If some ifl(x-l,y-j) fall out of image boundaries, replace them by 0 

//Do not consider reference pixels (I(x-i,y-j) where ij=0) 

//Apply (6) to compute e(x,y) 

sum_kl=0; 

for(i=0; i<a; i+=I) 

{ 

for 0=0;j<a;j+=l) 
{ 

if(I(x-hy'j)==out of bounds pixel) l(x-i,y'j)=0; 
if(i==0 && J==0) continue; 
sum_kl=sum_kl ^kg *I(x-i,y'j) ; 

} 

} 

sumjd =sum_ld> > 1 0; 
e(x,y) =I(x,y)-sum_kI; 

} 

} 



FIGURE 2 



t 



int RTJPEG::CompressRT(BYTE* input. BYTE* output, double ratio, double maxtime. 
int bpp=l, bool rgb=false) 

{ 

int sizejorig = Size(input); II original image size in bytes 
// Set max and min quality values 
intq0=5; // 5% quality 

int ql =95; II 95% quality 

// Find respective compression ratios. Higher ratio corresponds to lower quality 
double rO = sizejorig/Compress(input,output,qO, int bpp=l, bool rgb^false) ; 
II highest ratio 

double rl = size_orig/Compress(input,output,ql, int bpp=l, bool rgb^false); 
II lowest ratio 

// Return if proposed ratio is smaller than the smallest possible 
if(ratio<rl) return JPEG:: Compress(input, output, ql); 
II Return if proposed ratio is larger than the largest possible 
if(ratio>rO) return Compressfinput, output, qO); 
II Start timer 
clock J start, finish; 
start = clockQ; 

II Start iterative process to estimate the quality value 
double duration=0; double r - ratio; 

while (duration<maxtime) 

{ 

II Use linear quality estimate 

int q = qO+(r-rO)''(ql-qO)/(rl-rO); 

II Alternatively, a bi-section estimate can be used in the previous line: 
//intq = (qO+ql)/2; 

// Update estimated corresponding compression ratio value 

r = size_orig/Compress(input,output,q); 

II Update compression and quaUty boundaries 

if(ratio>r) 

{ 



} 

else 

{ 



rl=r; if(ql==<]) 6reaA:; // convergence 

= q: 



rO=r; if(qO==q) break; II convergence 

qO=q: 

} 

II Update timer 

finish = clockO; duration = (double) (finish - start) / CLOCKS J'ERJEC; 

} // end of iterative process 
// Compress 

return Compress (input, output, q); 

} 
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// It is assumed that ApplicationEntityList ael and DICOMObject d variables 
// have been aheady initialized based on the present network 
// configuration and the data to be sent. 

// Find current remote application entity, to which data must be sent 
ApplicationEntity ae = aelGetCurrentEntityQ; 

II Create PDU connection 
PDUpdu; 

II Start network performance timer 
NetworkTimer nt; 

II Find data size to be sent 

ntGetDataSize(d); 

nt.StartTimerQ 

II Send data to current application entity 
pduSend(d,ae); 

II Stop the timer, evaluate observed network bandwidth, and store it into ae 
nt. EstimateCurrentBandwidth(ae); 
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// It is assxuned that ApplicationEntityList ael and DICOMObject d variables 
// have been akeady initialized based on the present network 
// configuration and the data to be sent. 

// Compress DICOM object to accommodate current network bandwidth 
int dsize = d.GetSizeQ; 

int ntime = ae^GetDownloadTimeQ; II maximum allowable download time 

int ctime = dsize / ae.GetCurrentBandwidthQ; //time we'll spend with current bandwidth 

iffctime > ntime) II low bandwidth, we need compression 

{ 

if(ae. CanAcceptLossyO) II we can use lossy compression 
{ 

// Compress with maximum ratio allowed 

double ratio = min(ctime/ntime, ae.GetMaximumCompressionRatioQ); 
d. Compress(false, ratio); 

} 

else if(ae. CanAcceptCompressedQ) II only lossless, try our best 
{ 

d. Compress (true, /''^ignored*/ 1.0); 

} 

else II no JPEG compression allowed, do nothing 

{ 
} 

} 

else II network is fast enough, do not have to compress 

0 

II Create PDU connection 
PDUpdu; 

II Send compressed data to current application entity 
pdu.Send(d,ae); 
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class AccurateTimer 

{ 

private : 

int Initialized: 

int64 Frequency: 

int64 BeginTime: 

public : AccurateTimerQ // constructor 

{ 

II get the frequency of the counter 

Initialized - QueryPerformanceFrequency( (LARGE JNTEGER 
"^J&Frequency ): 

} 

void StartTimerQ II start timing 

{ 

if( ! Initialized ) return 0: II error - couldn't get frequency 
// get the starting counter value 

QueryPerformanceCounter( (LARGE JNTEGER V&BeginTime ): 
return 

} 

double EndTimerQ II stop timing and get elapsed time in seconds 
{ 

if( ! Initialized ) return 0.0: II error - couldn't get frequency 
// get the ending counter value 
int64 endtime: 

QueryPerformanceCounter( (LARGE JNTEGER "^J&endtime ): 
II determine the elapsed counts 

int64 elapsed = endtime - BeginTime: 

II convert counts to time in seconds and return it 
return (double) elapsed / (doubleJFrequency: 

} 
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bool SoundRecorder::Record(int sf bool stereo, int max_time=300, int 
max_size= 1000000, int format=WAV) 

{ 

DWORD dwRetum; 

II Open a waveform-audio device with a new file for recording. 
\i(!OpenMCI0) return; 

MCI_RECORD_PARMS mciRecordParms; 

II Set recording parameters here as fields in mciRecordParms 

mciRecordParms,dwFrom = 0; 

mciRecordParms. dwTo = maxjime; 

mciRecordParms.dwCallback = (DWORD) (this->GetSafeHwnd()); 
II Record 

if (dwRetum = mciSendCommand(mJDevice.wDeviceID, MCI_RECORD, 
MCI_FROM\MCI_NOTIFY, (DWORD) (LP VOID) &mciRecordParms)) 

{ 

AfaMessageBox(ErrorMCI(dwRetum, "Recording error: ")); 
mciSendCommand(mJDevice.wDeviceID, MCIjCLOSE, 0, NULL); 
return false; 

} 

return true; 

} 
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void SoundRecorder::Insert(DICOMObject& dob) 

{ 

BYTE* sound; 

if(GetFormatO != MPS) II we have to encode the sound 
{ 

SoundEncoder se; 

if(se.Encode(GetSoundQ, GetFormatQ, sound) 
{ 

return false; II we failed to encode 

} 

} 

//Store sound bytes into odb object 



return true; 

} 
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1 - llneamp Interpolation 

2 - parabolic root Inferpolation 
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